Correlation induced resonances in transport through coupled quantum dots 
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We investigate the effect of local electron correlations on transport through parallel quantum 
dots. The linear conductance as a function of gate voltage is strongly affected by the interplay of 
the interaction U and quantum interference. We find a pair of novel correlation induced resonances 
separated by an energy scale that depends exponentially on U . The effect is robust against a small 
detuning of the dot energy levels and occurs for arbitrary generic tunnel couplings. It should be 
observable in experiments on the basis of presently existing double-dot setups. 



The theoretical and experimental research on elec- 
tronic transport through ultrasmall quantum dots has 
become a very active field. Various fundamental physi- 
cal phenomena such as quantum interference, Coulomb 
blockade, and the Kondo effect, 0, 0] strongly affect 
the transport properties. Currently, the focus is shift- 
ing towards few-electron double-dot structures, that are 
studied as artificial molecules, interferometers, and for 
charge- and spin-based quantum computing, yj la, L3, killl 
flCM . HH Ha , Theoretical investigations into the role of 
correlations in systems of two or more coupled quantum 
dots are still at the beginning and much remains to be 
explored. In parallel quantum dots connected to common 
leads new physics is to be expected due to the interplay 
of correlations and quantum interference, an issue that is 
of interest also from a broader perspective. 

In the present Letter, we investigate this problem for 
a specific model of two parallel quantum dots coupled 
by an electron interaction U (see Fig. The effect 
of external electrostatic potentials on two-path inter- 
ference was studied earlier (magneto-electric Aharonov- 
Bohm effect). Q We here investigate the role of inter- 
action induced potentials. This is especially relevant for 
molecular transport, as the interaction induced poten- 
tials will be far larger than the external potentials which 
are difficult to apply to such small structures. We study 
how the linear conductance G as a function of gate volt- 
age V g changes with increasing U. Considering the entire 
parameter space we find a very rich generic behavior and 
predict the appearance of novel correlation induced reso- 
nances (CIRs) if U is larger than a critical interaction U c . 
The effect is robust: It appears for almost arbitrary com- 
binations of the four tunnel couplings and also remains 
visible for a small detuning of the dot level energies. The 
separation of the resonances in gate voltage defines an 
energy scale that depends exponentially on U and on a 
combination of the tunnel couplings. It is argued that 
this new correlation effect is unrelated to Kondo physics. 
We employ a powerful new method, the functional renor- 
malization group (fRG), to efficiently obtain both numer- 
ical and analytical results for this many-body problem, 
and we have confirmed all the essential features of our re- 
sults using the numerical renormalization group (NRG). 





FIG. 1: System of coupled quantum dots with common leads. 



Double-dot geometries that could form the basis to ver- 
ify our predictions have been experimentally realized in 
Refs. H B H Our model is equally of relevance for 
transport through two nearly degenerate levels of a single 
dot, a subject that has attracted much attention recently 
in attempts to understand the puzzling behavior of the 
transmission phase 0. 

We study a system of two quantum dots j — 1,2 each 
having a single level Sj as sketched in Fig. ^ The dots 
are coupled by a Coulomb interaction U > and are con- 
nected to two common leads I = L,R via tunnel barriers 
t l y The dot Hamiltonian is f/dot = J2j £ jd\dj + U{n\ — 
1/2) (ri2 — 1/2) and the dot-lead coupling is given by 
Ht = — J2j i(t l j c o i^j + H.c), where cj ; denotes the cre- 
ation operator at the end of the semi- infinite lead I. The 



leads are modeled by Hi = -tELoti/^+L' + H - c -)- 
The energy scale of the dot level broadening is given by 
= n\t l j\ 2 pi, where pi denotes the local density of states 
at the end of lead I. As usual we later take pi to be en- 
ergy independent (wide band limit). The ring structure 
is pierced by a magnetic flux <p that we take into account 



a i<t> 



2tt corre- 



by multiplying t% by a phase factor 
sponds to the flux quantum hc/e). For symmetry reasons 
one only has to consider < <fi < ir. The levels are shifted 
by a common gate voltage V g such that S\ = V g + 6 and 
t2 = V g , where S denotes a detuning. For 6 = due 
to the shift of rij in H^ ot , V g = corresponds to half 
filled dots. We neglect the spin of the electron and thus 
suppress the spin Kondo effect. Experimentally, the con- 
tribution of spin physics may be excluded by applying a 
sufficiently strong magnetic field. 01 We focus on tem- 
perature T = 0. The spectral properties 15] and level 
occupancies ^| of this model were investigated ear- 
lier. 
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To compute G and, in addition, the level occupancies the frequency-dependent conductance 
(rij) we mainly use a recently developed fRG scheme. p^j 



The starting point is an exact hierarchy of differential 
flow equations for the real-space self-energy matrix E A 
and higher order vertex functions, where A G (oo, 0] de- 
notes an infrared energy cutoff which is the flow param- 
eter. We truncate the hierarchy by neglecting the flow 
of the two-particle vertex only considering S A , which is 
then energy independent. This approximation and vari- 
ants of it were successfully used to study a variety of 
transport problems through quasi one-dimensional wires 
of correlated electrons (Tomonaga-Luttinger liquids) . 
In addition, spectral and transport properties of loc ally 
correlated systems were investigated using the fRG.|2fJ 
Although the guiding principle behind the above approxi- 
mation scheme is perturbation theory, it was shown that 
depending on the problem studied, the fRG procedure 
correctly leads to power-laws with U dependent expo- 
nents (Tomonaga-Luttinger liquids) as wen as expo- 
nential dependencies on U (Kondo effect). 



ReGH = -^|/ ,/| 2 <5(£/-£o-^), c^>0. (3) 



The current operator / = e(N R - N L )/2 (with Ni = 
v/ = i[H,Ni)) is expressed in terms 



Em=o c L,i c ™,i and N i 



of lead- and dot-operators. Its matrix elements Iqj (with 
indicating the ground state and / the excited states; 
Eq and Ef are the respective energies) are evaluated in 
the NRG basis and the ui — > limit then provides the 
conductance (7.|2ll| 

We first consider the case of degenerate levels with 
S = and later discuss how the results change for 5 > 0. 
The behavior of G(V g , U) can be cast in four classes that 
can already be identified at U = 0. (i) If two or more of 
the r^- are such that no closed path between the left 
and right lead exists G{V gi U) = 0. (ii) IyTf = r$T% 
and 4> — t: In this case one can introduce new fermionic 
dot states such that one only couples to the left lead and 



For the double-dot X A is a 2 x 2 matrix in the dot the other only to the right lead, implying G(V g , U) = 0. 

We note that in this case (and only in this case) a con- 
served pseudo-spin variable (left/right) exists. [13 (hi) A 



label j. The diagonal parts are real and V A = + V g 
can be considered as effective dot level positions. The 
(for arbitrary complex off-diagonal contribution t A = 
— Hi 2 is a hopping between the two dot states generated 
by the interaction. The flow equations are 
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with j being the complement of j. The Green function 
is Q A {iuo) = \ioj — h A (iuj)j with 

h A (iuj) : 



Vi A - iTi sgn(w) -t A - 17 sgn(w) 
- (t A )* - iY sgn( W ) V 2 A - iT 2 sgn( W ) 



(2) 



oi<t> 



and r, = x;,rj, 7 = vrTrf 

tial conditions are V A ~°° = V g + S, V{ 



f . The ini- 

= V g , and 

t A=ca = 0. To obtain an approximation for S and the 
one-particle Green function (using the Dyson equation) 
one has to solve the system Q of four real coupled dif- 
ferential equations. This can easily be done numerically 
and for a specific class of parameters also analytically. 
The occupancies (rij ) can directly be calculated from the 
Green function. For the present problem the same holds 
for G as current vertex corrections vanish. One can easily 
derive a lengthy expression for G in terms of the parame- 
ters r^- , (j> and the renormalized level positions Vj = V A=0 
and hopping td = t A=0 not presented here. Vj and td de- 
pend on T l j,4> as well as on V g , 5, and U. 

We use the NRG as a nonperturbative method to con- 
firm the essential validity of the physics discovered within 
the fRG. We allow arbitrary dot-lead couplings by em- 
ploying the Kubo- formula for the current, which yields 



nonvanishing (but nongeneric) conductance is found for 
rfTf = rfrf ^ and (j) = 0. In this case G(V g , U = 0) 
is given by a Lorentzian centered around 0. The U > 
dependence of G(V g ,U) can most easily be studied in 
the exactly solvable case of equal r' .lS] It is character- 
ized by two Coulomb blockade peaks located at w ±U/2. 
(iv) For all other Tj , <j>, that is for generic parameters on 
which we focus in the following, the peak in G(V g , U = 0) 
at V g = [as in case (hi)] is replaced by a dip with 
G(V g ,U = 0) = 0. For equal T l j, <f> = 0, S > 0, and 
U = the appearance of such dips was discussed earlier 
and explained as a destructive interference between path 
traversing dots 1 and 2 respectively. [2^ In the limit of 
a strong asymmetry of the transmission probability via 
dots 1 and 2, e.g. for T\ <C T l 2 , the dip can be viewed 
as a Fano anti-resonance resulting from the interference 
of a resonant path and a path with energy independent 
transmission: For energies at which the transmission via 
dot 1 shows a resonance, the transmission via dot 2 can 
be regarded as constant. 

Fig.[3shows the generic evolution of G(V g ) for increas- 
ing U at 6 — 0. Because of particle-hole symmetry G is 
symmetric around V g — 0. Energies are given in units 
of r = J2j Tj ■ Increasing U the height of the two peaks 
resulting from the dip at V g — increases and the maxi- 
mum flattens. At a critical U — U c ({Tj},4>) each of the 
peaks splits into two. For the present example the fRG 
approximation is U c /T w 4.69. Further increasing U the 
two outer most peaks move towards larger \V g \ and be- 
come the Coulomb blockade peaks located at V g « ±U/2. 
The other two peaks at ±Vcm are the novel CIRs, where 
Vcir > decreases with increasing U. Associated with 
G(V g = 0) = at U — is a jump of the transmis- 
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;Nll/r=4.69=U c /r- 




FIG. 2: (Color online) Generic results for G(V g )/(e 2 /h) (solid 
lines), (m) (dashed lines), and (7x2) (dashed-dotted lines) at 
different U obtained from the fRG with if = 0.27V, rf = 
0.33r, rf = 0.16r, = 0.24r, <f> = tt, and 5 = 0. The two 
novel correlation induced resonances are visible in the lower 
panels (large U), near V g = 0. 



sion phase by n. As the phase evolves continuously with 
U and particle-hole symmetry holds for any U, the ir- 
phase-jump and thus G(V g = 0) = must remain for all 
U. For U > U c the height of all four peaks is equal to 
h max ({Tj} , (f>) < e 2 /h and does not change with U. 

For equal T\, 4> = and 5 = [case (ii) above] the 
conserved pseudo-spin leads to orbital Kondo physics. (lij 
Remnants of this effect were found in the vicinity of this 
parameter point but die out quickly away from it. As our 
correlation effect appears generically, in particular for pa- 
rameters far away from the Kondo point it is apparently 
unrelated to Kondo physics. 

In addition to G(V g ) in Fig. [21 we present the dot oc- 
cupancies. For small U the V g dependence of the oc- 
cupancies of dots 1 and 2 is monotonic. In the oppo- 
site limit of large U, (rij) depends non-monotonically on 
Vg (HJE3 Starting at negative V g the level that is cou- 
pled more strongly (here j = 1) is depopulated at the 
first Coulomb blockade peak, while the occupancy of the 
other level stays close to 1. Close to V g = we find an 
inversion of the population followed by another depopu- 
lation of the more strongly coupled dot across the second 
Coulomb blockade peak at V g » U/2. Note that the 
non-monotonic behavior of the {rij) sets in for interac- 
tions smaller than U c and is thus not directly related to 
the appearance of the CIRs. This is consistent with the 
observation that in contrast to the CIRs (see below) the 
non-monotonicity of (rij) can already be observed within 
a self-consistent Hartree approximation. [To] 

The left part of Fig. |21 shows a comparison of fRG and 
NRG data for G at V g > and for different U. To clearly 
resolve the CIRs we use a logarithmic scale. The NRG 
data show all the features discussed in connection with 
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FIG. 3: (Color online) Comparison of fRG (lines) and NRG 
(symbols) results for the same parameters as in Fig. [5] Left: 
G(V g ) for different U (U = 3.5F: solid line and circles; U = 
7T: dashed line and squares; U = 14T: dashed-dotted line 
and diamonds). Right: The resonance position Vera as a 
function of U. 




FIG. 4: (Color online) Renormalized dot level positions Vi/T 
(dashed line) and V2/I (dashed-dotted line). The thick solid 
line is the product ViV^/F 2 . The thin solid line lies at 
— 4rfrf/T 2 and its intersections with the thick solid line 
determine the peak positions of G. The two crossings in- 
dicated by the arrows are at V g = ±Vcra- The parameters 
are U = 1QT, rf = rf 
5 = (U c /T « 5.05 from fRG) 



0.3r, I2 = I2 = 0.2r, = 7T, and 



Fig. |2l In particular, for increasing U > U c , Vein be- 
comes small quickly. The right part of Fig. shows the 
U dependence of Vcir extracted from the numerical data. 
For U sufficiently larger than U c , 



Vcm/T ex exp [-C ({v)} , 0) U/T] 



(4) 



with C > 0. By construction the fRG based approxi- 
mation scheme works particularly well for small to inter- 
mediate U. At larger U and for all Tpcj) we tested the 
fRG overestimates the positions of the Coulomb blockade 
peaks and the CIRs. As will be shown in an upcoming 
publication this can systematically be improved using a 
more elaborate fRG truncation scheme. For a specific 
class of Tj , <fi we next analytically confirm the exponential 
dependence of Vcir, on U and derive an explicit expres- 
sion for C using the fRG. In the most general case the 
dependence of C (and U c ) on r'- is complex and requires 
further investigation. Roughly speaking C increases (U c 
decreases) with increasing asymmetry of the [see also 
Eq. JSJ) below]. For fixed and increasing < <fi < tt, 
C decreases while U c increases. 

We now consider rf = rf, r£ = rf, but rf ^ rf, 
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FIG. 5: (Color online) Dependence of G(V g ) on the bare level 
splitting 8. The parameters are U = 3I\ rf = Q.7T, rf = 



o.2r, r£ = 

from fRG). 



o.02r, r 



R 



0.08r, and c/> = ([/ C /T » 0.842 



and <j> = 7T. In this case 7 = in Eq. (J2J and the off- 
diagonal elements of the Green function Q A are propor- 
tional to t^. Initially (at A = 00) t§ vanishes and it will 
thus remain zero during the fRG flow which leads to a 
simplification of the flow equations QJ. For small V^, 
that is small V gi these equations can be solved analyti- 
cally and in the limit U S> |Tf — T% | we obtain 



Vcm/r oc exp 



2tt 



1 2 



(5) 



In Fig. 0] we show the renormalized level positions Vj as a 
function of V g for fixed U > U c . For \V g \ > C7/2, Vi and 
V2 become equal and are given by V g — sgn(V^) U/2. For 
— U/2 ^ V g < the position of the level that is coupled 
more weakly (here j = 2) is smaller than the chemical 
potential /i = 0, while the other level has energy larger 
than [i. For < V g J$ U/2 the role of the two levels is 
interchanged. This explains the observed V g dependence 
of the (rij) discussed in connection with Fig. Peaks in 
G are found at V1V2 = — 4rf T% ■ In particular, the cross- 
ings of ViV 2 (thick solid line) and -4rfT2 (thin solid 
line) at ±Vcir occur because V\ and V2 continuously go 
through at V g = 0, one coming from above the other 
from below. This is not the case if the self-consistent 
Hartree-Fock approximation is used to compute the ef- 
fective level positions. Within this approach V1V2 discon- 
tinuously jumps at V g = leading to only two crossings 
at V g w ±U/2 associated with the Coulomb blockade 
peaks. This shows that similar to the Kondo effect more 
sophisticated methods than the Hartree-Fock approxima- 
tion are required to describe local correlation effects. The 
fRG is such a method (see also Ref. Efj) . 

We now investigate the stability of the CIRs in the 
presence of an initial level splitting S > 0. In Fig. 
we show fRG data for G(V g ) with different 6 at small 
\V g \. For S > T and all parameter sets we studied G only 
shows the two Coulomb blockade peaks at V g w —U/2 — S 
and V g k U/2. The way how this large 5 limit is reached 
depends on the specific choice of parameters but common 
to all cases, remnants of the CIRs were clearly observable 
for S <C r, as exemplified in Fig. |3J We found a similar 
S dependence using NRG. A direct hopping between the 



two dots induces a level splitting if the noninteracting 
dot Hamiltonian is diagonalized and has thus an effect 
similar to that of 5 > 0. 

To summarize, we found that the transport properties 
of a double-dot interferometer are strongly affected by 
electron correlations. We discovered a novel pair of cor- 
relation induced resonances that should be measurable 
in double-dots of appropriate geometry [tI IsL IsL [Toj in the 
presence of strong Zeeman splitting. Varying the tunnel 
barriers (and thus U/T) it should be possible to study 
the entire scenario discussed above. Apparently this cor- 
relation effect is unrelated to both spinQ and orbit alflij 
Kondo physics. Rather, it follows from the interplay of 
local correlations and quantum interference. It is thus 
likely that similar effects will appear in transport through 
more complex systems as e.g. ring-like molecules studied 
in the context of molecular electronics. Besides revealing 
interesting new physics, we showed that the fRG method 
is a very promising tool to investigate problems with local 
Coulomb correlations. In comparison to NRG the fRG 
is far superior in terms of the numerical effort required, 
e.g. enabling efficient analysis of parameter dependen- 
cies. Furthermore, the fRG can easily be extended to 
more complex systems with local electron correlations. 

We thank J. v. Delft, T. Enss, W. Hofstetter, Th. Pr- 
uschke, H. Schoeller, K. Schonhammer, A. Sedeki, and 
A. Weichselbaum for valuable discussions and in par- 
ticular M. Sindel for supplying us with his NRG code. 
V.M. is grateful to the Deutsche Forschungsgemeinschaft 
(SFB 602) for support. 



[1] 
[2] 

[3] 
[4] 
[5] 

[6] 

[7] 

[8] 

[9] 
[10] 

[11] 
[12] 

[13] 
[14] 



R. Schuster el at, Nature 385, 417 (1997); M. Avinun- 
Khalish et ai, Nature 436, 529 (2005), and Refs. therein 
Mesoscopic Electron Transport, edited by L.L. Sohn 
L.P. Kouwenhoven, and G. Schon (Kluwer, Dodrecht 
1997). 

L. Glazman and M. Raikh, JETP Lett. 47, 452 (1988) 
T. Ng and P. Lee, Phys. Rev. Lett. 61, 1768 (1988). 
D. Goldhaber-Gordon et ai, Nature 391, 156 (1998) 
W. van der Wiel et at, Science 289, 2105 (2000). 
D. Loss and D.P. DiVincenzo, Phys. Rev. A 
57, 120 (1998); D. Loss and E.V. Sukhorukov, 
Phys. Rev. Lett. 84, 1035 (2000). 

F. Marquardt and C. Bruder, Phys. Rev. B 68, 195305 
(2003). 

A. W. Holleitner et ai, Phys. Rev. Lett. 87, 256802 
(2001); Science 297, 70 (2002). 

J.C. Chen, A.M. Chang, and M.R. Melloch, 

Phys. Rev. Lett. 92, 176801 (2004). 

M. Sigrist et al, Phys. Rev. Lett. 93, 066802 (2004). 

A. W. Holleitner et ai, Phys. Rev. B 70, 075204 (2004). 

J.R. Petta et al, Phys. Rev. Lett. 93, 186802 (2004). 

F.H.L. Koppens et ai, Science 309, 1346 (2005). 

J.R. Petta et ai, Science 309, 2180 (2005). 

Yu.V. Nazarov, Phys. Rev. B 47, 2768 (1993); A. van 

Oudenaarden, M.H. Devoret, Yu.V. Nazarov, and 



5 



J.E. Mooij, Nature 391, 768 (1998); W.G. van der Wiel 

et al, Phys. Rev. B 67, 033307 (2003). 
[15] D. Boese, W. Hofstetter, and H. Schoeller, Phys. Rev. B 

64, 125309 (2001). 
[16] M. Sindel, A. Silva, Y. Oreg, and J. von Delft, 

Phys. Rev. B 72, 125316 (2005). 
[17] J. Konig and Y. Gefen, Phys. Rev. B 71, 201308(R) 

(2005). 

[18] V. Meden, W. Metzner, U. Schollwock, and K. Schon- 

hammer, J. Low Temp. Phys. 126, 1147 (2002). 
[19] V. Meden et al, Europhys. Lett. 64, 769 (2003); T. Enss 



et al, Phys. Rev. B 71, 155401 (2005); X. Barnabe- 
Theriault, A. Sedeki, V. Meden, and K. Schonhammer, 
Phys. Rev. Lett. 94, 136405 (2005). 

[20] R. Hedden, V. Meden, Th. Pruschke, and K. Schonham- 
mer, J. Phys.: Condens. Matter 16, 5279 (200 4); S. An - 
dergassen, T. Enss, and V. Meden, cond-mat/0509576 

[21] M. Sindel, W. Hofstetter, J. von Delft, and M. Kinder- 
mann, Phys. Rev. Lett. 94, 196602 (2005). 

[22] A. Silva, Y. Oreg, and Y. Gefen, Phys. Rev. B 66, 195316 
(2002). 



